Ab-initio calculation of the electronic and optical excitations in polythiophene: 

effects of intra- and interchain screening 
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We present an ab-initio calculation of the electronic and optical excitations of an isolated poly- 
thiophene chain as well as of bulk polythiophene. We use the GW approximation for the electronic 
self-energy and include excitonic effects by solving the electron-hole Bethe-Salpeter equation. The 
inclusion of interchain screening in the case of bulk polythiophene drastically reduces both the 
quasi-particle band gap and the exciton binding energies, but the optical gap is hardly affected. 
This finding is relevant for conjugated polymers in general. 
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Introduction - The interest in the electronic and op- 
tical properties of semiconducting conjugated polymers 
has increased enormously since the discovery of electrolu- 
minescence Q of these materials. The microscopic mod- 
elling of charge carriers and excitations in conjugated 
polymers is a subject of intensive research at present. 
In particular, knowledge of the single-particle bands and 
the electron- hole (exciton) binding energies is important. 
In conventional semiconductors such as GaAs the opti- 
cal excitations are well described in terms of very weakly 
bound electron- hole pairs (Wannier excitons) . In crystals 
made of small organic molecules such as anthracene, the 
electron- hole correlation energy is large, which essentially 
confines the exciton to a single molecule (Frenkel exci- 
ton). Exactly where conjugated polymers fit in between 
conventional semiconductors on the one hand and molec- 
ular crystals on the other, is the subject of a heated de- 
bate 1^ . Most of the experimental results have been ob- 
tained from optical experiments on the condensed (thin 
film or bulk) material. Most of the theoretical modelling 
has been at the level of a single molecule or a single poly- 
mer chain. In this paper we will focus on the differences 
between an isolated polymer chain and a bulk polymer 
material and show that these have a large effect on basic 
electronic parameters such as the band gap and the exci- 
ton binding energies. We calculate the electronic struc- 
ture and the optical excitation of conjugated polymers in 
an ab-initio manner, i.e. without 'tunable' parameters. 

Ab-initio calculations of several polymers within the 
Local Density Approximation of Density Functional The- 
ory (DFT-LDA) yield band gaps that are typically 40% 
smaller than the experimental optical gaps How- 
ever, it is well known that DFT-LDA formally cannot 
be applied to excited states. Moreover, no excitonic ef- 
fects are taken into account in these calculations. An ab- 



initio many-body GW approximation (GWA) Q calcula- 
tion on polyacetylene (PA) was performed by Ethridge et 
al. 1^ , but excitonic effects were not taken into account in 
their work cither. Very recently, a calculation by Rohlf- 
ing and Louie |^] for isolated chains of both PA and poly- 
phenylene-vinylene (PP V) , including both single- (quasi- 
particle) and two-particle (exciton) excitations, yielded 
optical gaps in good agreement with experiment. How- 
ever, the calculated lowest-lying singlet exciton binding 
energy for PPV in that work is much larger than in a 
recent experiment Q. 

In this article, we calculate single- and two-particle ex- 
citations for an isolated polythiophene (PT) chain as well 
as for bulk PT, where we will show that the differences 
between the two cases have very important consequences. 
We focus on PT, since it is one of the simplest and most 
studied polymers. Our main conclusions, however, are 
relevant for conjugated polymers in general. 

Computational Methods - Many successful ab-initio 
calculations of the quasi-particle (QP) band structure of 
conventional anorganic semiconductors have been per- 
formed with the GWA for the electronic self-energy S 
of the one-particle Green function. Very recently 
important progress has been made in the evaluation of 
the two-particle Green function, from which the optical 
properties can be obtained. This is done by solving the 
Bcthe-Salpeter equation ||lo| , |rT|] (BSE), which is equiva- 
lent to a two-body Schrodinger equation, for an electron 
and a hole forming an exciton. We adopt these two meth- 
ods and neglect the differences in geometry between the 
excited states and the ground state. DFT-LDA calcula- 
tions on oligothiophenes show that relaxation energies, 
calculated by relaxing an excited state from the ground 
state geometry to its optimized geometry, typically are in 
the range 0.1-0.2 eV Such values are upperbounds 
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for the corresponding states of the polymer. 

We start our calculations with a pseudo-potential 
plane-wave DFT-LDA calculation |^ of a geometry- 
relaxed FT chain in a tetragonal supercell (14.8 a.u. in 
the chain direction, 15 a.u. in the perpendicular direc- 
tions, the latter value large enough to consider the chains 
as isolated in this calculation). We use Hartree atomic 
units unless specified otherwise. The single-particle ex- 
citation energies Enk are evaluated by solving the QP 
equation 
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0nfc(r) 



J [Fpp(r,r') + E(r,r',S„fc)] (/)„fc(r')dV' = E„kMr), (1) 

where Vpp is the pseudo-potential of the atomic core, 
Vh the Hartree potential due to the valence electrons, 
and E the self-energy. We make the usual approx- 
imation that the QP wave functions </)„fc can be re- 
placed by the DFT-LDA wave functions. In DFT-LDA, 
I](r, r'jw) — 14c(r)(5(r — r'), with T4c(r) the exchange- 
correlation potential, while in the GWA the first term 
in the many-body expansion of S in terms of the Green 
function G and screened interaction W of the system 
is used, I](r,r',i) — iG{r,r',t)W{r,r',t), after Fourier 
transformation from frequency to time. The screened 
interaction is evaluated in the random-phase approxima- 
tion (RPA). We calculate these approximations of E and 
W in the mixed-space imaginary-time formalism |l3| , |l4| , 
from the DFT-LDA wave functions and energies. To 
study a truly isolated chain, we remove the 'crosstalk' 
between periodic images of the chain by taking as a unit 
cell the spatial region closer to the atomic positions of a 
specific chain than to those of any other; the Coulomb 
interaction v{r — r') = l/|r — r'| is then cut off by setting 
it zero if r and r' belong to different cells. 

The two-body electron-hole Schrodinger equation re- 
lated to the BSE is solved by expanding the cxci- 
ton wave functions <l'(re,r;j) of zero-momentum exci- 
tons, which are the optically relevant ones, in prod- 
ucts of conduction (c) and valence (v) wave functions, 
$(re,r,i) = I]fcct, ^fe™0cfe(re)</>*fc(r/O. The expansion co- 
efficients Akcv and the exciton binding energy Eb should 
obey 



[Eck — Eyk — Eg + Eb]Akc 



cv.k' c'v' 



KCV.K' c'v' 



0, 



(2) 



where Eg is the QP band gap, Eck and E^k are the QP 
energies, and Wkcv,k'c'v' and V^^^j,,^,^, (only present for 
singlet excitons) are the matrix elements of the static 
(uj = 0) screened interaction and the exchange matrix 
elements of the bare Coulomb interaction, respectively. 
Formally, dynamical screening effects may only be ig- 
nored in the BSE if Eg :$> E^. However, it has been 



shown that dynamical effects in the electron-hole screen- 
ing and in the one-particle Green function largely cancel 
each other and hence the static approximation is ex- 
pected to be accurate, even when the relation Eg ^ E^ 
does not strictly hold. 

As was shown by several authors [ p^p7| ], in a quasi- ID 
system (such as an isolated polymer chain with only in- 
trachain screening) there is no long-range screening. For 
realistic bulk polymers, however, both intra- and inter- 
chain screening are important, the latter giving rise to 
the screening at long distances. While it is in principle 
possible to perform a GM^A and exciton calculation for a 
crystalline structure, the amount of computational work 
is as yet prohibitively large. Instead, we choose to include 
the interchain screening in the following way, capturing 
the essential physics. From the imaginary- frequency de- 
pendent RPA polarizability tensor of an isolated chain 
in the long-wavelength limit we obtain a dielectric con- 
stant parallel (x) to the chain direction, e\\{iijj), and per- 
pendicular {y,z) to the chain direction, £_\_{iuj), by per- 
forming a 2D line-dipole lattice sum in the perpendicular 
directions, for an experimentally determined PT crystal 
structure | |l8| , p^ . Details of this approach will be given 
in a future article ||2^. We now take for the interchain 
screening interaction W^'^'^ = W — v the following form: 

{ [el{tLo)x' + e\\{tLo)e^{iu;){y' + z^)] ''^^ - v{r)} , (3) 

which is correct for distances much larger than the typ- 
ical interchain distance rjntor, for which we take 10 a.u. 
(typical for the crystal structure of Refs. |p^ , ^ ). The 
prefactor in this equation takes care of a soft cutoff for 
distances below renter , for which the interchain screening 
should vanish. We remark that the details of this interac- 
tion cannot be important since PT samples prepared in 
many different ways show very similar optical behavior; 
probably, most samples consist of small crystallites sep- 
arated by disordered regions. We now construct a total 
screened interaction 

Wtotai(r, r', iuj) = Wt^i,,{r, r', iuj) + 

Wi-'„(r-r',zc.) + i,(r-r'), (4) 

where Win"j.a ^he screening interaction already calcu- 
lated for the isolated chain. The interaction Eq. (jj) is 
correct at short and long range and is expected to be 
sufficiently good at intermediate range. We remark that 
the overlap between wave functions on, and therefore the 
electronic coupling between, neighboring chains is very 
small; hence, we can again use the isolated chain wave 
functions to calculate the Green function and self-energy. 
In our exciton calculations, we take the electron and hole 
on the same chain (so-called intrachain excitons). Hence, 
the only, but essential, difference between our calcula- 
tions for the isolated PT chain and bulk PT is in the 
screened interaction used. 
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We determine the various cutoff parameters in our cal- 
culations (number of /c-points, bands, grid mesh in space 
and imaginary time) such that the QP band gap is con- 
verged to within 0.05 eV. In the construction of the exci- 
ton wave functions we use only the tt and tt* wave func- 
tions, resulting in an accuracy of the lowest-lying exciton 
binding energies of about 0.1 eV. Below we give all cal- 
culated energies in eV with a precision of two decimal 
places. 

Results - The calculated GWA QP band structure of 
the isolated chain is shown in Fig. ^ We find a band 
gap of 3.59 eV (DFT-LDA: 1.22 eV). The binding ener- 
gies of the lowest-lying excitons and the resulting optical 
band gap are listed in Table |. The lowest-lying singlet 
exciton (^Bu) for the isolated chain has a binding en- 
ergy of 1.85 eV, leading to an optical gap of 1.74 eV. 
To our knowledge, no direct experimental information is 
available for either the QP band gap or the ^Bu exciton 
binding energy of PT separately. Only the difference, i.e. 
the optical gap, corresponding to the onset of the optical 
absorption in a well-ordered alkyl-substitutcd polythio- 
phene, is known to be about 1.8 eV 22|. This is in 
good agreement with our calculated optical gap. Simi- 
lar conclusions regarding the optical gap were reached by 
Rohlfing and Louie ^ for isolated PA and PPV chains. 

However, the difference between the ^Bu and ^ Ag bind- 
ing energies of the isolated PT chain is definitely not in 
agreement with a recent experiment p3], see Table |. 
Moreover, the ""^Bu exciton binding energy of 1.85 eV 
is very large compared to values currently discussed in 
the literature. There is a heated debate going on about 
the magnitudes of exciton binding energies in conjugated 
polymers Negligibly small (0.1 eV or les s p^ ), inter- 
mediate (~ 0.5 eV 1^), and large (~ 1.0 eV ]24[]) energies 
have been proposed. However, these values concern ei- 
ther films or bulk systems, both of which are essentially 
3D. 

The QP band structure calculated with our model 3D 
interaction Eq. (^) for bulk PT is also given in Fig. |^. 
The band gap has decreased to 2.49 eV. The results for 
the exciton binding energies are included in Table ^. The 
^Bu exciton binding energy has decreased to 0.76 eV, 
resulting in an optical gap of 1.73 eV. So, even though 
both the exciton binding energy and the band gap have 
changed considerably, the optical gap has hardly changed 
compared to the isolated chain. Furthermore, the rela- 
tive exciton energies are now in good agreement with the 
experimental data . If we fit our exciton coefficients 
to the hydrogcn-hkc form Ak/Ak=o ~ 1/(1 + Oex^^)^j 
we obtain for the ^Bu exciton sizes Ccx ^ 12 a.u. and 
^ 18 a.u. in the case of the isolated chain and the bulk 
situation, respectively. 

In order to test the sensitivity of our results to the 
precise value of the cutoff distance rintor we performed 
similar calculations for nnter = 8 a.u. and nnter = 12 a.u. 
The QP band gaps are 2.32 and 2.61 eV, respectively. 



The ^Bu binding energies are 0.64 and 0.86 eV and hence 
the optical gaps are 1.68 and 1.73 eV, respectively. This 
means that the optical gap is quite insensitive to the 
choice of rintcr- The energy differences between the exci- 
tons are also hardly influenced. 

Discussion and conclusions - Summing up, we have 
calculated the single-particle band structure and lowest- 
lying exciton binding energies of an isolated polythio- 
phene chain and bulk polythiophene. For the isolated 
chain (only intrachain screening) we find a large band gap 
and large exciton binding energies, due to the absence of 
long-range screening. Upon including interchain screen- 
ing, responsible for the long-range screening present in 
bulk polythiophene, we find that both the band gap and 
exciton binding energies are drastically reduced. How- 
ever, the optical gap hardly changes. 

We suggest that these conclusions hold for conjugated 
polymers in general. In this light, we can understand 
the fact that the calculations by Rohlfing and Louie 
on isolated chains of polyacetylene and PPV yield good 
results for the optical gaps, whereas their lowest-lying 
singlet exciton binding energy of 0.9 eV for PPV is far in 
excess of a recent experimental value of 0.35 ±0.15 eV 0, 
obtained by a direct STM measurement. The inclusion 
of interchain screening effects will drastically reduce the 
calculated binding energy (by more than a factor of two 
in our case of polythiophene), and may well lead to better 
agreement with this experiment. A similar drastic reduc- 
tion of the exciton binding energy by interchain screening 
effects was predicted recently by Moore and Yaron p^ ] 
in polyacetylene, within the semi-empirical Pariser-Parr- 
Pople theory. Clearly, it would be very interesting to 
repeat the experiment in Rcf. for polythiophene and 
polyacetylene. 

We conclude that a correct many-body description of 
the electronic properties of bulk polymer systems should 
include the effect of interchain screening. An impor- 
tant consequence is that neither Hartree-Fock nor DFT- 
LDA calculations should be relied upon in this context, 
since Hartree-Fock does not contain screening effects at 
all and since the exchange-correlation potential in DFT- 
LDA only depends on the local density and cannot de- 
scribe the non-local effects due to the long-range screen- 
ing. 

Financial support from NCF (Dutch National Com- 
puting Facilities) through project SC-496 is acknowl- 
edged. G. Brocks acknowledges the financial support 
from Philips Research through the FOM-LZM (funda- 
menteel onderzoek der materie - laboratorium zonder 
muren) program. We thank Michael Rohlfing for very 
fruitful discussions. 



3 



[1] J.H. Burroughes, D.D.C. Bradley, A.R. Brown, R.N. 

Marks, K. Mackay, R.H. Friend, P.L. Bum, A.B. Holmes, 

Nature 347, 359 (1990). 
[2] J.-L. Bredas, J. Cornil, A.J. Heeger, Adv. Mat. 8, 447 

(1996) . 

[3] G. Brocks, RJ. Kelly, R. Car, Synth. Met. 55-57, 4243 
(1993). 

[4] L. Hedin, Rhys. Rev. 139, A796 (1965); for a review of 
the GWA see: F. Aryasetiawan and O. Gunnarsson, Rep. 
Prog. Rhys. 61, 237 (1998). 

[5] E.G. Ethridge, J.L. Fry, M. Zaider, Rhys. Rev. B 53, 
3662 (1996). 

[6] M. Rohlfing and S.G. Louie, Rhys. Rev. Lett. 82, 1959 
(1999). 

[7] S.F. Alvarado, R.F. Seidler, D.G. Lidzey, D.D.C. Bradley, 

Rhys. Rev. Lett. 81, 1082 (1998). 
[8] S. Albrecht, L. Reining, R. Del Sole and G. Onida, Rhys. 

Rev. Lett. 80, 4510 (1998). 
[9] M. Rohlfing and S.G. Louie, Rhys. Rev. Lett. 81, 2312 

(1998). 

[10] L. Sham and T.M. Rice, Rhys. Rev. 144, 708 (1965). 

[11] G. Strinati, Rhys. Rev. B 29, 5718 (1984). 

[12] The properties of excited states in finite systems such as 
oligothiophenes can be calculated within the DFT-LDA 
formalism, provided they have a symmetry different from 
that of the ground state. For instance, the relaxation en- 
ergy associated with the ^Bu state of 12T (an oligothio- 
phene consisting of 12 thiophene rings) is 0.2 eV. Singlet 
states typically have a smaller relaxation energy. 

[13] H.N. Rojas, R.W. Godby, and R.J. Needs, Rhys. Rev. 
Lett. 74, 1827 (1995). 

[14] X. Blase, A. Rubio, S.G. Louie, M.L. Cohen, Rhys. Rev. 
B 52, R2225 (1995). 

[15] F. Bechstedt, K. Tenelsen, B. Adolph, R. Del Sole, Rhys. 
Rev. Lett. 78, 1528 (1997). 

[16] H.J. de Groot, P.A. Bobbert, and W. van Haeringen, 
Rhys. Rev. B 52, 11000 (1995). 

[17] H.J. Schulz, Rhys. Rev. Lett. 71, 1864 (1993) and refer- 
ences therein. 

[18] Z. Mo, K.-B. Lee, Y.B. Moon, M. Kobayashi, A.J. 

Heeger, F. Wudl, Macromol. 18, 1972 (1985). 
[19] S. Briickner and W. Rorzio, Makromol. Chem. 189, 961 

(1988). 

[20] J.-W. van der Horst, P.A. Bobbert, M.A.J. Michels, G. 

Brocks, R.J. Kelly, in preparation. 
[21] R. D. McCuUough, R. D. Lowe, M. Jayaraman, and D. 

L. Anderson, J. Org. Chem. 58, 904 (1993). 
[22] K. Sakurai, H. Tachibana, N. Shiga, C. Terakura, M. 

Matsumoto and Y. Tokura, Rhys. Rev. B 56, 9552 

(1997) . 

[23] T.W. Hagler, K. Pakbaz, A.J. Heeger, Rhys. Rev. B 49, 
10968 (1994). 

[24] M. Chandross, S. Mazumdar, S. Jeglinski, X. Wei, Z.V. 
Vardeny, E.W. Kwock, T.M. Miller, Rhys. Rev. B 50, 
14702 (1994). 

[25] E.E. Moore and D. Yaron, J. Chem. Rhys. 109, 6147 

(1998) . 




-25 I 1 

r k X 



FIG. 1. The single-particle band structure of the valence 
electrons of polythiophene along the chain direction calcu- 
lated within the GWA, including intrachain screening only 
(isolated chain, dashed) and including both intra- and inter- 
chain screening (bulk, full lines). The DFT-LDA highest va- 
lence (tt) and lowest conduction (tt*) bands are drawn dotted. 
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TABLE L Binding and transition energies of the four low- 
est-lying excitons in polythiophene, and quasi-particle {Eg) 
and optical {Eo) band gaps calculated with intrachain screen- 
ing only (isolated chain) and intra- plus interchain screening 
(bulk). Experimental data before/after lattice relaxation ac- 
cording to ReL lH). All in eV. 
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